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Abstract 

To integer programming problems, computational algebraic approaches using Grobner bases or 
standard pairs via the discreteness of toric ideals have been studied in recent years. Although these 
approaches have not given improved time complexity bound compared with existing methods for 
solving integer programming problems, these give algebraic analysis of their structures. In this 
paper, we focus on the case that the coefficient matrix is unimodular, especially on the primal 
and dual minimum cost flow problems, whose structure is rather well-known, but new structures 
can be revealed by our approach. We study the Grobner bases and standard pairs for unimodular 
programming, and give the maximum number of dual feasible bases in terms of the volume of 
polytopes. And for the minimum cost flow problems, we characterize reduced Grobner bases in 
terms of graphs, and give bounds for the number of dual (resp. primal) feasible bases of the primal 
(resp. dual) problems: for the primal problems the minimum and the maximum are shown to be 
1 and the Catalan number 2( 2 ^d~P)> wnne f° r the dual problems the lower bound is shown to 
be Sl(2 L d / 6 J ) . To analyze arithmetic degrees, we use two approaches: one is the relation between 
reduced Grobner bases and standard pairs, where the corresponding relation on the minimum cost 
flow — between a subset of circuits and dual feasible bases — has not been so clear, the other is 
the results in combinatorics related with toric ideals. 

1 Introduction 

Recently, some algebraic approaches to integer programming problems have been studied. The two 
main approaches are using Grobner bases || and standard pair decompositions ||. Although they 
neither give improved complexity bounds compared with existing methods nor have been demonstrated 
to solve hard practical instances which cannot be handled by existing methods, these approaches 
themselves are very interesting by applying computational algebraic methods to such hard problems, 
and give algebraic analysis of structure of integer programming problems ||, [8L [| 17, 21, 22, p3[ . 



For an ideal over a polynomial ring, the reduced Grobner basis and the set of standard pairs are dual 
in a sense that the complement of the monomials in the initial ideal, which is generated by initial 
terms of the reduced Grobner basis, is the set of standard monomials, whose nice decomposition is 
the standard pair decomposition. This kind of duality may shed new light on duality in combinatorial 
optimization, and by considering a nice subclass of integer programming problems where the duality 
theorem holds, we might be able to obtain some complexity bounds by making use of the characteristics 
of the subclass, which could not be derived for general integer programming problems. 

The problems whose coefficient matrices are unimodular form a nice subclass in a sense that the 
system yA < c becomes totally dual integral (TDI). Then each standard pair corresponds to a dual 
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feasible basis, and the method using standard pairs is equivalent to calculate the reduced cost for each 
basis (Theorem |3.4[ ). Thus, the number of standard pairs, which is equal to that of dual feasible bases, 
gives the complexity of this approach. Additionally, the maximum number of standard pairs can be 
described by the normalized volume of another matrix (Theorem 

Especially, the minimum cost flow problems form a well-known subclass of unimodular integer 
programming problems which can be solved in polynomial time. A Grobner basis approach for the 
minimum cost flow problems is a variant of cycle-canceling algorithm. In the case of the strongly 
polynomial time algorithms || [l2|, |l3|| , for any feasible flow they choose polynomial size of negative- 
cost cycles (by the selecting rules) from the set of negative-cost cycles in the residual network, which 
may be of exponential size, as many as possible. Similarly, the algorithm using Grobner basis calculates 
the optimal flow by augmenting flows along the negative-cost cycles which correspond to the elements 
of Grobner basis. Thus the cardinalities of reduced Grobner bases may give some time bound for 
this algorithm. On the other hand, a standard pair approach for the minimum cost flow problems 
first finds the set of standard pairs, and solves linear system of equations for each standard pair until 
an integer and non-negative solution is obtained. For a network optimization problem, the duality 
between the reduced Grobner basis and the set of standard pairs corresponds to the relation between 
circuits and dual feasible co-trees, dually, cutsets and primal feasible trees. Since such a relation has 
not been so clear, the computational algebraic duality may be interesting method for the analysis of 
network problems. 

This paper is organized as follows. In Section 2, reduced Grobner bases and standard pairs 
are defined, and their relations with integer programming problems, regular triangulations and dual 
polyhedra are introduced. The case that the coefficient matrix is unimodular is studied in Section 3. 
The maximum arithmetic degree (i.e. the maximum number of dual feasible bases) is shown to be the 
normalized volume of a polytope defined by homogenizing the coefficient matrix (Theorem |3.3| ). In 
Section 4, we study the Grobner bases and standard pairs on the primal minimum cost flow problems 
on acyclic tournament graphs with d vertices. We show that three types of reduced Grobner bases can 
be characterized in terms of the circuits (Theorem 4J3, L£, 4.1C| ). These examples give the minimum 
and the maximum number of dual feasible bases of the minimum cost flow problems: the minimum is 



1 (Theorem 4.14) and the maximum is g( 2 ^_i ) (Theorem 4.15| ). This maximum is shown using the 



result in Section 3 and the result about the hyper geometric system on unipotent matrices and related 
polytope ||. In Section 5, we study the dual minimum cost flow problems. One reduced Grobner 
basis is characterized in terms of cutsets (Theorem |T^). We also show that the lower bound for the 
number of primal feasible bases of the minimum cost flow problems is f2(2L rf / 6 J). 





Grobner basis [§] 


Standard pair [£] 


Primal 


Term on graph 


Set of circuits 


Set of spanning trees 


Algorithm 


Variant of 
cycle-canceling 


Enumeration of 
dual feasible bases 


On acyclic tournament 
graph with d vertices 


min : d(d - l)/2 
max : ? 


min : 1 

max: M 2i tP) 


Dual 


Term on graph 


Set of cutsets 


Set of co-trees 


Algorithm 


Variant of 
cutset-canceling 


Enumeration of 
primal feasible bases 


On acyclic tournament 
graph with d vertices 


min : d — 1 

max : ? 


Lower bound 0(2L rf / 6 J ) 



Table 1: Dual algebraic approaches for primal and dual minimum cost flow problems 
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2 Toric ideals and Grobner bases 



For a matrix A G Z dxn and a cost vector c € M n , let IPa,c be the family of integer programming 
problems IPa jC (6) := minimize {c ■ x\Ax = b, x G N 71 } as b varies in {.Am|m G N 71 } C (N is the 
set of non-negative integers). The cost vector c is called generic if each program in IPa,c has the 
unique optimal solution. 

Let A; be a field and k[x] := k[xi, . . . , x n ] the polynomial ring. For an exponent vector a = 
(ai,...,a n ) G N n , we denote A total order on monomials in A; [a;] is a torm order 

if 1 is the unique minimal element, and >- x v implies x u+w y x v+w for all u,v,w G N n . For 
a fixed term order >-, the refinement >-c of c by >- is a total order such that tc' u >-c x v if either 
c • u > c • v or "c • w = c • i> and x u >- x v " holds. If c > 0, then >~c becomes a term order. 

The toric ideal I a of A is a binomial ideal I a '■= {x u — x v \ Au = Av, u,v G N n ). For any 
/ G Ia, the initial term iriy c (f) of / is the largest term in / with respect to >~c- Then we define the 
initial ideal inyJJ.A) of 1a as iny c {lA) '■= {i n )~ c {f) I / € Ia)- 

2.1 Grobner bases and Conti- Traverse- algorithm 

A finite subset Qy c = {51, • • • ,g s } Q Ia is a Grobner basis for Ia with respect to y c if iny c (lA) = 
(iriy c (gi), . . . ,iny c (g s )}. In addition, Grobner basis Qy c is reduced if £/^ c satisfies that (i) for any i, 
the coefficient of iny c (gi) is 1, and (ii) for any i, any term of is not divisible by iny c (gj) (i 7^ j). 
If is a term order, then the reduced Grobner basis Qy c exists uniquely, and is calculated by 
Buchberger algorithm (see Q). Any Grobner basis for I a is a basis of I a J1- 

Ia is called homogeneous with respect to the positive grading deg(xj) = di > (i = 1, . . . , n) if, 
for any / = /1 + /2 + ■ • • + f m G (/j is the homogeneous component of degree £ in /), /j G I4 for 
any £. Then I4 is homogeneous if and only if I a is generated by homogeneous polynomials 0. 



Proposition 2.1 ([|16|) If Ia is a homogeneous with respect to some positive grading deg(xj) = di > 
0, then >~c becomes a term order for any c G W 1 \ {0}, and the reduced Grobner basis Gy c exists. 

The support supp(it) of a vector u is the index set {i \ Ui 7^ 0}. Any u G TP" can be written 
uniquely as u — u~^~ — u where , it G BJ n and have disjoint support. Then Gyp can be written as 
Qy c = {x u i — x u i , . . . , x u p — x u p } for some finite u\, . . . , u p G ker(yl) n Z n Q. 

Example 2.2 Ie£ A = ( ^ J and consider the minimum cost flow problem IPA,c(b) = 



-10 1, 

minimize{c ■ x \ Ax = b,x = {x\^,x\ 3,^2,3) G N 3 }. T/ien t/ie toric idea/ is Ia = (^1,2^2,3 — ^13)- 



2 O O 3 

Figure 1: Acyclic tournament graph with 3 vertices. 



If c = (01,2,01,3,02,3) = (3,1,2), then inc{lA) = (^1,2^2,3} and reduced Grobner basis is Qy c = 
{^1,2^2,3 - ^1,3}- 

In the rest of this paper, we consider a cost vector c which >-p becomes a term order for some term 
order >-. Let IPA,y c (b) be the problem to find the unique minimal element in {x G N n [ Ax = b} with 
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respect to y c . Then the solution u of IP^y^ip) is one of the optimal solutions of IPA,c(b). Conti 
and Traverso M introduced an algorithm based on a Grobner basis to solve IPA,y c (b). We describe 
the condensed version of the Conti- Traverso Algorithm (see |l6| ) , which is useful for highlighting the 
main computational step involved. The normal form of / € k[x] by the reduced Grobner basis Q is 
the unique remainder obtained upon dividing / by Q. 

Algorithm 2.3 (Conti- Traverso Algorithm |2], |l6| ]) 

1. Compute the reduced Grobner basis Qy c of I a with respect to >-c- 

2. For any feasible solution v of IPa,c(P), compute the normal form x u of x v by Gy c - 

3. Output u. u is the optimal solution of IPA,y c {b). 

Thus reduced Grobner basis Q)~ c is a test set for IPa^ c [|(], 21 1. 



Example 2.2 (continued) Let b = (4,5,-9). For a feasible solution (4,0,9), the normal form of 
x i,2 x 2,3 by Q>~ c is #i 3^2,3- Thus the optimal solution of IPA,y c {b) is (0,4,5). 



2.2 Standard pair decompositions 

Let [n] := {1, . . . , n}. For a monomial x a G k[x] and an index set a C [n], (x a , a) is a standard pair of 
iny c (lA) if (i) supp(a) n a = 0, (ii) every monomial in x a ■ k[xj \ j G a] := {x a ■ f \ f G k[xj \ j G a]} 
is not an element of m^ c (I/i), and (hi) there exists no other (x a ,cr'), which satisfies (i) and (ii), 
such that x a> divides x a and supp(a; a /a; a ') Ua C a'. We denote S(in^ c (lA)) the set of all standard 
pairs of in,y c (lA)- We use the same (x a ,a) to denote the set of all monomials in x a • k[xj \ j G a]. 
Then the above condition (hi) says that (x a ,o-) <fi {x a ,o-') for any other ) which satisfies the 

condition (i) and (ii). The standard pairs of m^ c (/^) induce a unique covering for the set of standard 
monomials of iriy c (lA), which we call the standard pair decomposition of m^ c (J^). LS'(in>- c (/ i 4)) is 
called the arithmetic degree of iriy c (lA) and denoted by arith-deg (iriy c (lA)) jl8| . 

Example 2.2 (continued) Forc= (3,1,2), the standard pairs of in (3,1,2) (1 a) ar e (1> {(1> 2), (1, 3)}) 
and (1, {(1, 3), (2, 3)}), thus the arithmetic degree of in (3,1 ,2) (1 a) is 2. On the other hand, for c = 
(1, 4, 2), the standard pair o/m( 1) 4, 2 )(/A) is (1, {(1, 2), (2, 3)}), thus the arithmetic degree of inn ,4^) (Ja) 
is 1. 




Figure 2: Two types of standard pair decompositions. A point (p, q, r) in the figure corresponds to 
the monomial x\ 2 X \ 3^2 3- 



Let c be a generic cost vector. Then iny c (lA) = Wc(-Ta)- Let {01, . . . , a n } be the column vectors 
of A and cone (A) the cone generated by a\,...,a n . For a C [n], we denote A a for the submatrix 
of A whose columns are indexed by a. For a cost vector c, we define the regular triangulation Ac 
of cone (A) as follows: cone(A CT ) is a face of Ac if and only if there exists a vector y G M. d such that 
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y ■ dj = Cj (j G a) and y ■ aj < Cj (j ^ a). If cone(A a ) is a face of A c , a also is called a /ace of A c . 
The genericity of c implies that Ac is in fact a triangulation (i.e. each face of Ac is simplicial) |[7]]. 



Lemma 2.4 (gg, gg) 

(i) If inc(lA) has (*,cr) as a standard pair, then a is a face of Ac- 

(ii) inc{lA) has (1,<t) as a standard pair if and only if a is a maximal face of Ac- 

(iii) If ai, ... ,a n span an affine hyperplane, then Ac 2s i/ie same as £/ie regular triangulation of 
conv(^4) with respect to c, and the number of standard pairs (*, a) for a maximal face a of Ac 
equals the normalized volume of a in Ac- 



When vertices of conv(A) are in the m-dimensional lattice L ~ Z m , we define the normalized 
volume of a maximal face a of Ac by the volume of a with the normalization that the volume of the 
convex hull of 0, ex, ... , e m is 1. Here, {ej}i<j< m are the basis of the lattice L. 

For a polyhedron P C R n and a face F of P, the normal cone of i 7 at P is the cone Np{F) := 
{lo € K n | uo ■ x' > u; • x for all x' E F and x G -P}. The set of normal cones for all faces of P is called 
the normal fan of P. 

Lemma 2.5 ([^]) Ac is £/ie normal fan of the polyhedron Pc '■= {y €M. d | yA < c}. 

We remark that Pc is the dual polyhedron for the linear relaxation problem LPA,c{b) := minimize 
{c • x | yla; = 6, cc > 0} of IPA t c(b). When A is row-full rank, this lemma shows that there is one-to- 
one correspondence between the dual feasible bases of LPA,c{b) and the maximal faces of Ac- 



Example 2.2 (continued) For c= (3,1,2), A 



(3,1,2) 



{{1,2},{2,3},{1},{2},{3},0}. 





Figure 3: Dual polyhedron ^(3,1,2) an d regular triangulation A( 3 



Let u be the optimal solution to IPa,c(P). Since standard pairs cover w&c(/a)j x u is covered 
by some standard pair (x a ,a). Thus u = a + Yliea ^i e i f°r some non-negative integers {/cj}j ecr , 



and b = Au = A (a + Ylt^a ^i e i) = Aa + Yliea ^i a i- Lemma 2A implies that {ai}i &a are linearly 
independent. Therefore {k.j}i ecr is the unique solution to the linear system ^ieo- x * a « = b — Aa. This 
observation induces an algorithm to solve IPA,c(b) using the standard pair decomposition of ztic(Za)- 

Algorithm 2.6 (Solving IPA,c{b) using S(inc(lA))) 

(i) For (x a ,o~) G S{inc{lA)), solve the linear system X^eo-^* * = b — Aa. Let {ki}i^ a be the solution. 

(ii) If {ki}i£ a are both integral and non-negative, output a + X^ieo- ^» e « as ^ e optimal solution. Oth- 

erwise, repeat (i) for another standard pair. 

This algorithm solves at most arith-deg (inc(l yi))-many linear systems of equations. Therefore 
arith-deg (otc(-Ta)) is a measure of the complexity of IPa,c- 
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3 Standard pairs for unimodular programming 

Let A G 7j dxn be row-full rank and unimodular, i.e. each non-zero maximal minor is ±k for some 
k G N. Then WcC^a) is minimally generated by square-free monomials for any c fl6|| , and all standard 
pairs are obtained from all maximal faces of Ac- 

Lemma 3.1 ([]9|]) Let {mi,... ,ra s } be the minimal generators of inc(lA)- If mi, . . . ,m s are all 
square-free then S{inc{lA)) = {(L ') I °~ is the maximal faces of A c }. 

For a matrix A G Z dxn , the homogenized matrix A' G z( d + 1 ) x ( n + 1 ) of A is 

a'-=( 1 1 ••• 1 m = t 1 1 ■" 1 m a) 

V A J \ ai a 2 ■■■ a n J' {L) 

Let a! i = f _ ) for 1 < i < n and be the (n + l)-th column vector of A'. We remark that 

a'i, . . . , a^, a' n+1 span an affine hyperplane. 

We define another family IPa',(c,o) °f integer programming problem 

I^' l( c,Q) (&,/?) :=m™ze {c ■ x | A'^J = Q , (^J G N" +1 } 

as (^) varies in {A'u \ u G N n+1 }. We remark that (c, 0) is generic if c is generic. 

The next proposition is due to Sturmfels et al. |L8| for general ideals. We give another proof for 
the case of toric ideal. 

Proposition 3.2 (0]§) (x a ,a) G S(in c (I A )) (x a G k[x\, a C [n]) i/and on/y i/ (a; a , a U {n + 1}) G 
S (in (c,o) (I A'))- 



Proof: We first show that any monomial in (x a , a) is standard for inc(lA) if and only if any monomial 
in (x a , a U {n + 1}) is standard for in( C Q) (-^4')- Suppose that any monomial in (x a , a) is standard for 
inc(lA) and choose any x u x^ l+l G (x a ,a U {n + 1}). If there exist any other (Y) G N n+1 such that 
A'd) = A'(™) and ( v t ) ^ then = ,4t>, and cu<cv since cc u ^ !n c (J A ). Therefore, (™) is 
the optimal solution to IPa' ,(c,o){Au, Y^h=i u i + ^ there does not exist such (Y), then clearly (^) 
is the optimal for IPa> ,(c,o)(^ u > S^=i M « + k). This shows that any monomial in (x a , a U {n + 1}) is 
standard for ^( Ci0 ) (//!')• 

Conversely, suppose that any monomial in (x a , a U {n + 1}) is standard for m( C0 ) (^A') and choose 
any a? 11 G (a; a ,cj). If there exists some v G N n such that At? = Aia, then -A'(^) = ^.'(g) f° r any 
non-negative integers (/ such that p — q = XT=i u « ~~ X^=i n *- Since £E u :e^ +1 G (a3 a , <r U {n + 1}), 
(c, 0) • ) < (c, 0) • ( v ) , which implies that c ■ u < c ■ v. Therefore, u is the optimal solution to 
IPa,c(Au). If there does not exist such v, then clearly u is the optimal for IPa,c{Au). Thus any 
monomial in (x a ,a) is standard for mc(-Ol)- 

Let (x a ,a U {n + 1}) G 5(in( Cj0 )(/ J 4'))- K (a^f) C (:r a ',T) for any other (x a ' ,r) which sat- 
isfies (i) and (ii) in the definition of standard pairs for inc(lA), then (x a , a U {n + 1}) must be 

contained in ^a; a , r U {n + 1}^ , which contradicts the assumption. Thus (x a , a) G S(inc(lA))- On 
the other hand, if (x a , a U {n + 1}) ^ S(in^ c ^(Ia')), then there exists some ^ai a x^ +1 , t'^J which 
contains (x a ,aU {n + 1}) and ^a; a x^ +1 ,r'^ satisfies (i) and (ii) in the definition of standard pairs 

for in( C ) (Ia 1 ) • Then n + 1 G r', and therefore k = 0. Therefore, [x a ' ,t' \ {n + 1}) contains (a3 a ,cj) 
and satisfies (i) and (ii) in the definition of standard pairs. Thus (x a , a) ^ S(inc{lA))- This completes 
the proof. □ 
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Example 2.2 (continued) For this A, enlarged matrix A' is 

/ 1 111 

a! = 1 ioo 
\ -l o i o 

We consider I^i C fc[a;i,2> ^1,3, #2,3) a^]. For c = (3,1,2), the standard pairs of W(3,i 2,o)(^A') are 
(1, {(1, 2), (1, 3), 4}) and (1, {(1, 3), (2, 3), 4}). On the other hand, for c = (1, A, 2), the standard pairs 
°/ ^(1,4,2,0) CTa') are (1) {(1) 2), (1, 3), (2, 3)}) and (1, {(1, 2), (2, 3), 4}). In this case, the only standard 
pair (1, {(1, 2), (2, 3), 4}) satisfies the condition in Proposition \3.tj which corresponds to the standard 
pair (1, {(1, 2), (2, 3)}) of in {1A2) {I A ). 

Since a^, . . . , a' n+1 span an affine hyperplane, the normalized volume of conv(^4') gives the number 
of standard pairs of in^c,k)(^A') which correspond to maximal faces of A(cfc) by Lemma |2.4j (hi). 

Theorem 3.3 ([|7|]) If A is a unimodular matrix, then the maximum arithmetic degree of inciJ-A) 
equals the normalized volume of conv(A'). 



Proof: For any c, the set of standard pairs oEi%c(Ia) is {(1 °~) | <7 is a maximal face ofAc}, and each 
(1,<t) corresponds to the standard pair (1, a U {n + 1}) of ini^Q\{I A i). Especially, a U {n + 1} is a 
maximal face of A', co y Therefore, 

arith-deg (in c (lA)) = |{(1 a) G S (in c (I A ))}\ 

= \{(l,o-U{n + l})eS(in ic>0) {lA>))}\ 

|(*,r) G S (w(c,0)(-^4')) I T : maximal face of A'^ c0 ^| 
normalized volume of conv(j4'). 



< 



Let Ia C k[x] and I A ' C k[x\, . . . ,x n ,x n+ i]. Then x a — x^x^ +1 G Ia> (x a , x^ G k[x\) if and 

only if Y17=i( a i ~ bi) = k and x a — x b G I a- We consider that c = (l,l,...,l) and >~ is any reverse 
lexicographic term order such that x n+ \ is the smallest variable. Then for any g in the reduced 
Grobner basis Q for I A > with respect to ^(c,o)j % (Co) (s) does not contain x n+ \ by the definition of 
the term order, and in^ (Co) {g) is square-free since {"V (Co) (<?) I 9 G minimally generates in^/ (1^4) 

for some term order y'. Thus the corresponding triangulation A(_ is unimodular [p^| , and each 

(c,o) 

facet of A(_ Cq corresponds to a standard pair of inc(lA) injectively. Then the arithmetic degree of 
incijA) is equal to the number of facets of A(_ , which is the normalized volume of conv(^f). □ 

(C,0) 

We consider the primal problem which is equivalent with LP^c(6) : 

P(M I) city := maximize {(— c) T a/ | Mx' + i^a;" = bs, x , a;" > 0}, 
which corresponds to some basis 73, and its dual problem 

^(7 -mt) := minimize i b Bv" I - M T y" = c, y',y" > 0}, 

where M G Z dx (™~ d ), 5 = (b B ,&*) (6b = (>i)feB G Z d , = (6*)^b = G Z n ~ rf ), J d G Z dxrf 
and I n -d G ;g( n_c x ( n_(i ) are identity matrices, x" (resp. x') is a basic (resp. non-basic) variable for 
P(M I) c(6)) y' (resp. y") is a basic (resp. non-basic) variable for D ^ ^(c), and c is a reduced 
cost for B. 
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For any standard pair (1, a) of m c (I^) = mg(/m n), a := {1, . . . , n} \ a forms a basis of the dual 
problem D ^ ^(c) (Lemma ^5| ). Let 

cri := ({1, . . . ,n} \ B) Her, a 2 :=Bf]a, a\ := ({1, . . . ,<n}\B) Da, := B Da. 

Then the reduced cost of D „ r (c) for the basis a is 

(/ -M T ),£r ' 

b a = b a - NjiB- 1 )^, where B 1 = (-M T )^), Nt = (I ai (-M T ) ff2 ). 



Theorem 3.4 The solution of the equation in Algorithm \2.(\ (i) for a standard pair (1,<t) is the 

reduced cost of D m ? Cc) for the basis W. 
J (I -M T ),b v ' J 

Proof: We show that b a is a solution of the linear system in Algorithm ^1] (i) for (1, a), i.e. (M CTl I a2 )b a = 
b. This is because 



(M ai I U2 )b a = (M ai I U2 )b a - (M ai I a2 )Nt{B^) L b<f 

= I a2 b a2 - (M ai (I ai ) T + / CT2 ((-M T ) (T2 ) T )( J B 1 - 1 )% 

= I a2 b a2 — {(M — A%fe) T ) + (-M - /^((-M T W) T )} (B^fb? (2) 

= I n b n + (A% %)(% (_M T )^) T ( J Br 1 ) T ^ 

= / CT2 6 CT2 + I w )Bj(B^fh 

= Ia 2 ba 2 + (MoT I-^)bw 

J(j 2 60-2 + I-^b^ 

= b, (3) 

the equation © follows from the fact that M = MI = M ai (I ai ) T + A%(I 5T ) T and —M = I(-M) = 
/ CT2 ((-M T ) (72 ) T + /^((-M T W) T . 
□ 



4 Grobner bases and standard pairs of the primal minimum cost 
flow problems 

Let Gd be the acyclic tournament graph with vertices 1,2, ... ,d and n = („) arcs, where each arc 
(i,j) (i < j) is directed from i to j. We consider the following minimum cost flow problem P^cfi): 

PA,c(b) '■= minimize {c T x \ Ax = b, x > 0}, 

where A € Z rfxra is the vertex-arc incidence matrix of Gd- 

A walk in Gd is a sequence (t>i,t>2, • • • ,v p ) of vertices such that (vi,Vi+i) or (vi+i,Vi) is an arc of 
Gd for each 1 < i < p. A cycle is a walk (v±, v%, . . . , v p , vi). A circuit is a cycle (ui,U2, ■ ■ ■ , v p , vi) such 
that Vi / Vj for any i ^ j- 

Definition 4.1 Lei C be a circuit in Gd and fix a direction of C. If C passes an arc (i,j) ufj times 
forwardly and iQ times backwardly, then we define Uq = (Ujj)i<i<j<d, Uq = (itZ)i<i<.;<d € t". The 
vector uc '■= u^, — Uq is called the incidence vector of C . We identify a cycle C of Gd with the 
binomial fc '■= x u c — x u c g I A . 



8 



Definition 4.2 A non-zero vector u in ker(A) is a circuit if its support supp(u) is minimal with 
respect to inclusion and the elements of u are relatively prime. When u £ ker(j4) is a circuit, we also 
call x u+ — x u a circuit of I a- We denote Ca ol set of circuits of I a- 

Then Ca corresponds to the set of all circuits in Gd- Let IAa be the union of all reduced Grobner 
bases for I a with respect to all term orders, which is called the universal Grobner basis of 1a- 

Proposition 4.3 ( ||16| ) For the vertex-arc incidence matrix A of Gd, Ua = Ca- Especially, any 
reduced Grobner basis of I a is square-free, and the number of elements in IAa is of exponential order 
with respect to d. 

Proposition 4.4 1a is not homogeneous for the grading deg(ajjj) = 1, but is homogeneous for the 
grading deg(xij) =j-i. 



Proof: For any d, x\p,xi,z — ^1,3 € 1a and £1,2^2,3 £ Ia- This implies that 1a is not homogeneous for 
the grading deg(xjj) = 1. 

Let vijt^j ■■■ ,v p ,vi be a circuit in Gd, C + := {k \ v k < v k+ i} and C~ := {k \ v k > v k+ i} (we set 
v p+1 := vi). The binomial f c corresponding to C is f c = ]lfcec+ x v k v k+1 - TLeC- x v k+1 v k - Then Ze- 
is homogeneous for the grading deg(xjj) = j — i because 

- deg Yl x v k +iv k = ( Vk + l ~ v k)- Yl ( Vk ~ Vk +^ 

\kec- J kec+ kec- 

v 

= ^(v k+ i - v k ) = 0. 




k=l 



□ 



Thus reduced Grobner basis exists for any c € R n \ {0} by Proposition 2.1 



4.1 Some Grobner bases for the primal problem 

We first show that the elements in reduced Grobner bases with respect to some specific term orders 
can be given in terms of graphs. As a corollary, we can show that there exist term orders for which 
reduced Grobner bases remain in polynomial order. For other applications of the Grobner bases found 
in this section, see our paper 1TT[ . 



Proposition 4.5 Let y be the purely lexicographic order induced by the variable ordering such that 
x i,j <~ x k,i if and only if 'i < k or (i = k and j < I). Then the reduced Grobner basis for I a with respect 
to y is {g ijk := x it jXj jk - x ijk \ % < j < k} U {g ijk i ■= x itk x jt i - x i} iXj )k | i < j < k < /}. In particular, 
the number of elements in this Grobner basis is equal to ( 3 ) + (^) . 

The set {gij k \ i < j < k} corresponds to all of the circuits of length three, and {gij k l \ i < j < k < 1} 
corresponds to some circuits of length four uniquely determined for each four vertices i, j, k, I (Figure ||). 



Proof: By Proposition |4.3| , it suffices to show that any binomial which corresponds to a circuit in Gd 
is gijk, dijki or whose initial term is divisible by some iny(gij k ) or iriy(gij k i). 
Any binomial corresponding to a circuit of length 3 is contained in {gijk}. 

The circuits defined by four vertices i < j < k < I are C\ := (i, j,k,l,i), C2 ■= (i,j,l,k,i), 
C3 := (i,k,j,l,i) and their opposites. The binomial which corresponds to C\ or its opposite is 
±(xijXj^ k x k j — Xij), whose initial term X{jXj k x k 1 is divisible by in^(g; L j k ). Similarly, the initial term 
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3 



Figure 4: Circuit corresponding to gij k (left) and circuit corresponding to g^jy (right). 



of binomial which corresponds to C2 or its opposite is divisible by iny^g^ji). The binomial which 
corresponds to C3 or its opposite is gijki- 

Let C be a circuit of length more than five. Let v\ be the vertex whose label is minimum in C, 
and C := U2, ■ ■ ■ , v p , Vi). Without loss of generality, we set V2 < v p . Let fc be the binomial 
corresponding to C, then iny(fc) is the product of all variables whose associated arcs have the same 
direction as (wi,^) on C. If V2 < V3, then (ui,«2) and (^2,^3) have the same direction on C. Thus 
both x„ lj1;2 and x V2tV3 appear in iriy(fc), and iny(fc) is divisible by in^(g 1 , 1 ^ 2 „ 3 ). If V2 > V3, then 
since v$ < V2 < v p , there exists some k (3 < k < p — 1) such that i>i < v k < V2 < Wfe+i- Then both 
x Vli „ a and x VktVk+1 appear in iriy(fc), and iriy(fc) is divisible by ^(fl^ 1^2^+1)- a 

Theorem 4.6 Zei >- be any term order and c = (ci^, . . . , ci^, 02,3, . . . , Cd-i.d) € R n satisfy (Hj+Cj^ > 
c i,k for any i < j < k and Ci tk + Cj i > Cij + Cj k for any i < j < k < I. Then the reduced Grobner basis 



for I a with respect to >-c is the same as the basis in Proposition l^.l 



Proof: Let >-' be the term order defined in Proposition ^4.5| . Then iny c (gij k ) = XijXj jk = in>y'(9ijk) 
since aj + c Jjfc > c i>k , and iny c (g ijk i) = x iyk x^i = iny>(g ijk i) since c iik + Cjj > + c j>k . Thus 
iriy c (lA) = iriyr(Iji), which implies that the reduced Grobner bases for I a with respect to >-c and >-' 
are the same. □ 

Proposition 4.7 Let >- be the purely lexicographic order induced by the variable ordering such that 
Xij >~ x k j if and only if j — i < I — k or (j — i = I — k and i < k). Then the reduced Grobner basis for 
I A with respect to y is {g ijk := x itj x jjk - x ijk \ i < j < k} U {g ijk i := Xijx j;k - x ijk Xjj \ i < j < k <l}. 
In particular, the number of elements in this Grobner basis is equal to ( 3 ) + (^) . 

The set {g%j k \ i < j < k) corresponds to all of the circuits of length three in Gd, and {gijki \ i < 
j < k < 1} corresponds to the set of circuits of length four in Figure || but the direction is opposite. 

Proof: Any binomial corresponds to a circuit of length 3 is contained in {gij k }- 

The circuits defined by four vertices i < j < k < I are C\ := (i,j,k,l,i), C2 ■= (i,j,l,k,i), 
C3 := (i,k,j,l,i) and their opposites. The binomial which corresponds to C\ or its opposite is 
±(xijXj jk x k j — Xij), whose initial term XijXj tk x k j is divisible by iriy(gij k ). The binomial which 
corresponds to C2 or its opposite is ±(xijXjj — Xi ik x k j). If its initial term is XijXjj, it is divisible by 
iny(giji). On the other hand, if initial term is X{ ik x k j, it is divisible by iny(gi k i). The binomial which 
corresponds to C3 or its opposite is g^u- 

Let C be a circuit of length more than five. Let (t>i,t>2) {v\ < V2) be the arc in C which the 
difference of labels is minimum, and C := (i>i,i>2, • • • ,v p ,vi). Let fc be the binomial corresponding 
to C, then iriy(fc) is the product of all variables whose associated arcs have the same direction with 
(vi,v 2 ) on C. 

If V2 < V3, then both x VltV2 and x V2tV3 appear in iny(fc), and iriy(fc) is divisible by iriy(g VlV2V3 ). 
Similarly, if v p < v±, then iriy(fc) is divisible by iny(g VpVlV2 ). 

Let U3 < V2 and v± < v p . Then v$ < v\ < V2 < v p by the definition of v\ and v%- If there exists 
some q such that v q < v q+ \ < v q+ 2, then iriy(fc) is divisible by iny(g VqVq+lV +2 ). Consider the case 



10 



that there does not exist such q. For any s such that v s < V\ < v s+ \ < v 2 , v s+2 < v s+ i by assumption, 
and v s + 2 < v\ by the definition of v± and v 2 . Thus there must be some r (3 < r < p — 1) such that 
v r < v\ < v 2 < v r+ i since v 3 < v 1 < v 2 < i p . Then in y (fc) is divisible by in y (g VrVlVaVr+1 ). □ 



Theorem 4.8 Let y be any term order and c = (ci j2 , • • • , c\ d, c 2,3> • • • , Cd—xd) £ ^ n satisfy Ci t j+Cj k > 
c i,k f or an y i < j < k and Cj^ + Cjf. > Ci k + Cj t i for any i < j < k < I. Then the reduced Grobner basis 
for I a with respect to is the same as the basis in Proposition |^, % 



Proof: Let y' be the term order defined in Proposition [4.7| . Then iny c {gij^) = Xi jXj^ = in y i{gijk) 
since + c jtk > Cj |fe , and in yc (g ijk i) = x^Xj^ = in y >(g ijk i) since c i){ + c jtk > c ijk + c,y. Thus 
m^ c (/A) = ot^'(/a)) which implies that the reduced Grobner bases for I a with respect to yc and y' 
are the same. □ 



Proposition 4.9 Let y be the purely lexicographic order induced by the variable ordering such that 
x i,j 'y x k,i if and only if i < k or (i = k and j > I). Then the reduced Grobner basis for La with respect 
to y is {gij := j | i < j — 1}. Ln particular, the number of elements in this 

Grobner basis is equal to (f) — (d — 1). 

The set {gij \ i < j — 1} corresponds to all of the fundamental circuits of Gd for the spanning tree 
T := {{i,i + 1) | 1 < i < d}. 

Proof: Let C be a circuit which is not a fundamental circuit of T. Let v\ be the vertex whose label 
is minimum in C, and C := (i^,^, • • • ,v p ,v±). Without loss of generality, we set v 2 < v p . Then the 
variable x VliV appears in the initial term of the associated binomial fc- Thus in y (fc) is divisible by 

in y (g VlVp ). □ 



Theorem 4.10 Let y be any term order and c = {c\ 2, ■ ■ ■ > c\ d> c i 3j 



,Cd-i,d) S M n satisfy c i:j > 



+ Cj-!j for any i < j 



Cj,i+1 + Q+l,«+2 + 

yc is the same as the basis in Proposition U71. 



1 . Then the reduced Grobner basis for La with respect to 



Proof: Let y' be the term order defined in Proposition [O]. Then in yc (gij) = Xij = in y i(gij) since 
Cij > c^,^+l+c^+l,^+2 + ■ ■ ' + c j-i,j- Thus in^ c (/yi) = in y i(LA), which implies that the reduced Grobner 
bases for La with respect to yc and y' are the same. □ 



4.2 Bounds for the size of Grobner bases 

Generally the degree of any reduced Grobner basis for toric ideal is of exponential order with respect 
to the number of rows in the matrix pq] , but the cardinality is not well understood. For the case of the 
toric ideals of acyclic tournament graphs, since those vertex-arc incidence matrices are unimodular, 
the cardinalities of the reduced Grobner bases may be bounded. 

Proposition 4.11 The minimum cardinality of the reduced Grobner bases for La is (2) — (d — 1). 

The basis we have shown in Proposition is the example achieving this cardinality. 



Proof: Since the reduced Grobner basis forms a basis for La, the cardinality of the reduced Grobner 
basis is more than that of the basis for La- Since La corresponds to the cycle space of Gd, the 
cardinality of the basis for La is equal to the dimension of the cycle space, which is (2) — (d — 1). □ 

To analyze the upper bound for the cardinalities of the reduced Grobner bases, we calculate all 
reduced Grobner bases for small d using TiGERS [10]. Table § is the result for d = 4, 5, 6, 7. 
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d 


# GB 


max cardinality 


min cardinality 


4 


10 


5 


3 


5 


211 


15 


6 


6 


48312 


37 


10 


7 


> 37665 


> 75 


15 



Table 2: Number of reduced Grobner bases, maximum and minimum of cardinality. 



For the case of d = 7, the number of reduced Grobner bases and the maximum of the cardinality 
are both too large, so we could not know the correct values. For d < 5, the reduced Grobner basis in 



Proposition 4.5 is the example achieving maximum cardinality, but for d > 6 the maximum cardinality 



is a little larger than the cardinality of Grobner basis in Proposition |4.5| . For d = 6, we do not know 
what cost vectors produce the Grobner bases of cardinality 37. The reduced Grobner bases which 
achieve the maximum cardinality seem to be complicated and difficult to characterize. 

Problem 4.12 Are the cardinalities of reduced Grobner bases for I a of polynomial order with respect 
to d? 

4.3 Standard pairs for primal problem 

In this section, we assume that c is generic. If c is not generic, then we may use d which is obtained 
by perturbing c such that d is generic and in c (f) contains the term in c >(f) for any / 6 I a- Since 
one constraint of PA,c(b) is redundant, we can consider the problem Pj c (b), which is obtained from 
Pa,c(P) by deleting the last constraint. Then mc(i^) = OTc(-f^)> and A is row-full rank. In addition, 
the regular triangulation of cone(^4) and that of cone(^4) by c are the same as a simplicial complex, 
thus we denote both triangulations A c . 

Since any initial ideal inc(lA) is generated by square-free monomials (Proposition |4.3| ), the set of 
standard pairs S{inc{lA)) are (1)°") where a ranges among all maximal faces of Ac- 

Let E be a set of arcs in Gd- For S C E, we denote x s := Y\(i j)es x i,j- 

The arcs in the optimum flow of uncapacitated minimum cost flow problem define a forest Q. 
Therefore, with the fact the dimension of cone(^4) equals d — 1, the next proposition is implied by 



Lemma 2.4, Proposition 4.3 and Lemma 3.1 



Proposition 4.13 (x a ,a) is a standard pair of inciJ-A) if and only if x a = 1 and a is a spanning 
tree of Gd such that x a ^ inc(lA)- 

Because of the result in Section |3|, there is one-to-one correspondence between the standard pairs 
(1, *) of inc{lA) and the dual feasible bases of P-^ c (b). Therefore, Algorithm ^1] for the minimum 
cost flow problem P^cip) is a variant of the enumeration of dual feasible bases. 

Grobner bases which have shown in the previous section give upper and lower bounds for the 
arithmetic degree (i.e. bounds for the number of vertices of the dual polyhedron). The genericity of 
c implies that the arithmetic degree of inc{lA) is equal to or greater than 1. 

Theorem 4.14 The minimum arithmetic degree of inc(lA) which c varies all generic cost vectors 
equals 1. 



Proof: For a cost vector c as in Theorem 4.10 , inc(lA) = (xij \ j — i > !)• Then x a £ inc(lA) 
if and only if a^j = for any (i,j) such that j — i > 1. The set of all such monomials coincides 
(1, {(1, 2), (2, 3), ...,(<£ — 1, d)}). Thus only this pair is a standard pair of inc(lA)- d 
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Theorem 4.15 The maximum arithmetic degree of iric{lA) which c varies all generic cost vectors 
equals Cd-i '■= ^C^d-i )' w hich is the (d — l)-th Catalan number. 

To show this theorem, we use the next result due to Gelfand et al. || which studies about some 
hypergeometric function. 

Lemma 4.16 ([j^]) Let A' be the enlarged matrix (||) for the incidence matrix A of the acyclic tour- 
nament graph with d vertices, and conv(^L') be the convex hull of a' 1: . . . ,a' n+1 . Then the normalized 
volume of conv(A') equals the (d — l)-th Catalan number Cd-i- 



Proof of Theorem 4-15 : Since A is unimodular, arith-deg (inc(lA)) < (normalized volume of conv(^4')) 



Cd-i by Theorem 



Because of Proposition 4.15 and Theorem [4.6| , for c as in Theorem |4.6| , (1, a) is a standard pair of 
inc{lA) if and only if a is a spanning tree of the acyclic tournament graph which satisfies the following 
two conditions: 

(a) there are nol<i<j<fc<d such that both (i, j) and (j, k) are arcs in a, and 

(b) there are no l<i<j<k<l<d such that both (i, k) and (j, /) are arcs in a. 

The number of such spanning trees are known to be the (d — l)-th Catalan number (e.g. see |14]]). □ 
We remark that the Catalan number equals ^j n 3/2 (l + O (-)) (e.g. see |J). 



5 Grobner bases and standard pairs of dual minimum cost flow 
problems 

In this section, we analyze Grobner bases and standard pairs for the dual minimum cost flow problems. 
As in Section ||, we study the problem which is equivalent with P^ jC (6): 

P, M ^ g(b) := maximize {(— c) T x' \ Mx' + Ix" = bs, x',x" > 0}, 

which corresponds to the basis {(1, 2), (2, 3), . . . , (d — 1, d)}, and its dual problem 

D _ AfT) ^(c) := minimize {b B y" \ Iy' - M T y" = c, y' , y" > 0}, 

where (M I) (resp. (/ — M T )) is the fundamental cutset (resp. fundamental circuit) matrix which 
corresponds to the spanning tree {(1_, 2), (2, 3) , . . . , (d—1, d)}, c isj:he reduced cost corresponding to the 
basis {(1, 2), (2, 3), . . . , (d - 1, d)}, b = (b B , b N ) = (&ij)i<j<j<d, b B = {bi,i+i)i<i<d, b N = (&ij)i<j-i = 
0, and 

X = (x', x"), x' = (Xi, 3 . . . , X ltd , X 2 ,4, • • • , Xd-2,d), x" = (xi >2 , X 2 ,3, X d -l,d) 

y = (y",y'), y' = (2/1,3 • • - ,2/1,^,2/2,4, • • -,yd-2,d), y" = (2/1,2,2/2,3,- • -,yd-i,d)- 

Then P (M /) 6 (6) has d - 1 constraints (i.e. (M I) G Z^ d -^ xn ), D _ A/T) ^(c) has n - d + 1 con- 
straints (i.e. (I - M T ) G z( n - d+1 ) xn ). 

Let Gd = (V, E). D C E is a cutset in Gd if there exists a partition (V\, V2) of F (i.e. V\ fl V2 = 
0, Vx UV 2 = V) such that D = {(i,j) e E \ i G V x and j G V 2 , or i G V 2 and j G Vi}. 

Definition 5.1 Let D be a cutset in Gd which corresponds to V = (V + ,V~). We define the vector 
u D G W 1 as 

{1 (i G V + and j G V~ ) 
-1 (i G V' and j G V+) . 
(otherwise) 

The vector up is called the incidence vector of D. 
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We identify a cutset D which corresponds to {V + ,V~) with the binomial fjj := x u d — x u d. 
Since the rank of the fundamental circuit matrix (/ — M T ) is n — d + 1 and each row vector of the 
fundamental cutset matrix (M I) is in ker((J -M T )), the set of row vectors of the fundamental cutset 
matrix (M I) forms a basis of ker(J — M T ). 

For the fundamental circuit matrix (I — M T ), the set of circuits C^j _m t ) corresponds to the set of 
all cutsets of G^- Since the fundamental circuit matrix (/ — M T ) is totally unimodular (e.g. see ]19|]), 
Proposition [P^ implies Crj _m t ) = Un -m t )- 

Proposition 5.2 For a cost vector b such that the linear system (M I)x = fog has a non-negative 
solution, I^j has a reduced Grobner basis with respect to b. 

Proof: Let a > be a solution of (M I)x = fag. We denote r» the i-ih. row of (M I), i.e. the row 
which corresponds to the fundamental cutset for the arc (i,i + 1). For each cutset D corresponds to 
(V+,V\ V+) (V+ C {1, . . . ,d - 1}), since u D = EieV+, i+i$V+ r i ~ E^y+, ;+ieV+ r ^ 

au D = ^ ari- ^2 a-Vi 

i&V+, i+l£V+ i£V+, i+l&V+ 

iev+, i+i(£v+ i<£v+, i+iev+ 

= b • u D . 

Thus m a (/o) = in?(fD) for any cutset D, and in a (I(i _m t )) = ^^(i — M T ) ) • Since a > °> -A/ T ) 
has a reduced Grobner basis with respect to b. □ 

Example 2.2 (continued) Lei c = (3,1,2) and fo = (4,5,-9). Then the primal and dual problem 
which corresponds to the spanning tree {(1,2), (2,3)} are the following. 

min 4yi )2 + 9y 2 ,3 

' J/1,3 

s.t. (11-1 -1)1 J/12 



s.i. 



1,3 






1 


1 o \ 


/ #1,3 


1 


o i ; 


a?l,2 






V x 2,3 



J/2,3 



»1,2, »1,3, «2,3 > J/l, 2 , 2/1,3, J/2,3 > 



T/ien Jq j -1,-1) = (a?i,2 — ^2,3, ^1,2^1,3 — 1, ^1,3^2,3 — 1) and reduced Grobner basis for b = (4,0,9) is 
{X2,3 ~ ^1,2, £1,2^1,3 - 1}. 

5.1 Grobner basis for dual problems 

As for primal problems, we show that the elements in reduced Grobner basis to some specific term 
order can be given in terms of graphs. 



Theorem 5.3 Let b be the cost vector which satisfies the condition in Proposition 5.1, fej.i+i > kjj+i 
(1 < v i < v j < d) and bij = ( v i, j such that j > i + 1). Then the reduced Grobner basis for 
I(i -A/T) with respect to b is { gi : — Y[j<i x j,i ~ llfc>i x i,k I * — 2,3,..., dj- . In particular, the number 
of elements in this Grobner basis is equal to d — 1 . 

Each gi is an incidence vector of the cutset which corresponds to (V \ {i}, {i})- 

Proof: For a cutset D which corresponds to (V + , V~) such that 1 E V + , we define P + := {i £ V + \ i 7^ 
d, i+ 1 G I/"} and := {i G | i ^ d, i+1 G Let P + = {h, . . . , i p } (n < i 2 < ■ ■ ■ < ip) 
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and P~ = {j'i, . . . , j q } (ji < j 2 < ■■■ < jq)- Then p = qorp = q + l, and i x < j x < h < h < * ■ • < 
*A < Jfc < ik+i < Jk+i < • • • Since 5 • w+ = 6 irjir+ i > £^ =1 = b ■ w r» ^(/d) = 

Since inr(g h+ i) = l\j< h Xj,h+i, m^/o) can be reduced by inAg il+1 ). □ 



5.2 Bounds for the size of Grobner bases 

Proposition 5.4 T/ie minimum cardinality of the reduced Grobner bases for In _j\/t) is d — 1. 77ie 
frasis we have shown in Theorem \5.3{ is the example achieving this cardinality. 



Proof: Since the reduced Grobner basis forms a basis for In _m t )> * ne cardinality of the reduced 
Grobner basis is more than that of the basis for In _ m t), which is d — 1. □ 

To analyze the upper bound for the cardinalities of the reduced Grobner bases, we calculate all 
reduced Grobner bases for small d using TiGERS [10]. Table S is the result for d = 4, 5, 6, 7. 



d 


#GB 


max cardinality 


min cardinality 


4 


7 


5 


3 


5 


48 


10 


4 


6 


820 


20 


5 


7 


44288 


39 


6 



Table 3: Number of reduced Grobner bases of dual problems, maximum and minimum of cardinality. 

We do not know what cost vectors produce the Grobner bases of maximum cardinality. The 
reduced Grobner bases which achieve the maximum cardinality seem to be complicated and difficult 
to characterize. 

Problem 5.5 Are the cardinalities of reduced Grobner bases for In _m t ) °f polynomial order with 
respect to d? 



5.3 Standard pairs for dual problem 

In this section, we assume that b is generic same as Section [4.3| . Since any initial ideal inriln _m t )) 
is generated by square-free monomials, any standard pair in S(inr(In _ m t^)) is of the form (1, *). 
Moreover, since the support of each optimal solution of D,j _ m t\ does not include a cutset, with 
the fact that dimcone((/ — M T )) = n — d + 1, the next proposition is implied by Lemma 2.4 and 



Lemma 3.1 



Proposition 5.6 (x a , a) is a standard pair of inriln _m t )) an d on ^V if xCL = 1 and o~ is a co-tree 
of Gd such that x a ^ inr(In _m t ))- 

Example 2.2 (continued) For c = (3, 1, 2) and b = (4, 5, —9), the initial ideal ^(4,0,9) (-^(ll— 1,-1)) = 
(^2,3, &i 2^1,3} has two standard pairs (1,{(1,2)}) and (1,{(1,3)}). 



Theorem 5.7 For any b which satisfies the condition in Proposition 5A, there exists S C {1, ... ,d — 
1} with \S\ > \_(d—l)/G\ such that, for any a C S, there exists a spanning tree T a of Gd which satisfies 
the following: 

(A) T a contains the arc set {(i,i + 1) j i G S\a} and does not contain any arc in + 1) | j £ a}, 
and 
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(B) (l,T a ) is a standard pair of inr(I^j _jv/ t ))> where T a := E \ T a is a co-tree ofT a . 

Especially, since T a ^ T T for any a, r C S (a ^ r), mr(I(j ~m t )) ^ as °^ least f7 (2 L rf / 6 J ) standard pairs 
for any generic b which satisfies the condition in Proposition \5.$\ . 



Proof: We divide {1, . . . , d — 1} into the following subsets. 



M 


:= {i 


€{1,.. 


.,d- 


1} 


Xij+i e m^(J (J _ AfT) )} 


Mi 


:= {i 


e{i,.. 


.,d- 


1} 


i i M , i = (mod 3)} 


M 2 


:={• 


€{1,.. 


.,d- 


1} 


i (£ M , i = 1 (mod 3)} 


Ms 


:={* 


€{1,.. 


.,d- 


1} 


i M , % = 2 (mod 3)} 



Lemma 5.8 |M | < f(d- l)/2]. 



Proof of Lemma 5.8: We consider a cutset L> which corresponds to (V + ,V ) such that fjj contains 



hi 



as a term of degree 1. Without loss of generality, we set i G V + . We assume that j — i > 1. Then 
for any A; (i < k < j), if k G ^ + then /d contains x^j and Xjj in the same term, otherwise fjj contains 
Xik and Xij in the same term, which contradicts that Xi.j is a term of fo of degree 1. Thus j = i + 1. 
In addition, G y~ for any A; < i and G y + for any k > i + l. Therefore, V + = {i, i + 2, i + 3, . . . , d} 
and V~ = {1, ... ,i- + 1}. 

We consider that inr(/) = aJi,»+i for some / G I(j _m t )- K G inr{I^ -m t ))i then / can be 

reduced by the binomial corresponding to the cutset between {i — 1, i + 1, . . . , d} and {1, . . . , i — 2, i} 
to 




;i+2 

and its initial term is Xa+i. Since both terms of this binomial contain Xjj+i, this implies that 



mr (/'/a^j-fi) = 1. Since b defines a term order by Proposition 5.2 , this is a contradiction. 

Similarly, Xj+i,i+2 £ inrily _m t ))- Thus |M | < f(d - l)/2]. □ 

Thus at least one of Mi, M2, M3 has at least \_{d — l)/6j elements. Let S be one of such 
Mi (i = 1,2,3). For any a := {i\ > 12 > • • • > i r } C S, we construct desired spanning trees 
T 0> T {ii}> T {n,; 2 }' T °- inductively. 

• Initial step: 

Let T$ := {(1, 2), (2, 3), . . . , (d — 1, d)}. Clearly T0 is a spanning tree. Since the reduced Grobner 
basis corresponds to a subset of cutsets, the initial term of any elements of reduced Grobner basis 
contains a variable x^j+i for some i. Thus x Tll> ^ mr(I(j _m t ))- 

• Induction step: 

Let T a \u\ be the desired spanning tree for a\{i r }. We define two edge set 

T 1 := {T ctUM \ {(i r , v + 1)}} U {(i r , i r + 2)} , T 2 := {T 1 \ {(i r + l,i r + 2)}} U {(i r -l,i r + 1)} . 

Then both T 1 and T 2 are spanning trees and satisfy the condition (A). We show that either T 1 or T 2 
satisfies the condition (B). 

(a) The case that T 1 satisfies the condition (B). Then T 1 is a desired spanning tree T a . 

(b) The case that T 1 does not satisfy the condition (B). 

In this case x T1 G inr(I^j -m t ))- Let Q be the reduced Grobner basis for _m t ) with respect to b. 
Then x Tl can be reduced some binomial g G Q, and such g is one of the following form (See Figure ||). 
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(i) <7?j which corresponds to the cutset for (V + , V ), V + = {p,p + 1, . . . , i r ,i r + 2, i r + 3, . . . , d} 

and V~ = {1, 2, . . . ,p — 1, i r + 1} for some p < i r , and its initial term is a product of variables 
corresponds to arcs from V + to V~ , or 

(ii) (The case of r > 1) <?L t ) which corresponds to the cutset for (V + ,V~), V~ = {1,2, . . . ,p — 

1, i r + 1, i q{l) + 1, . . . , i q (t) + 1} and V + = V \ V~ for 1 < 3 q(t) < ■■■ < 3 q(l) < r such that 
{iq{k) + 1) *<j(fc) + 2) G To-^i^} for k = 1, . . . ,t and 1 < 3 p < i r , and its initial term is a product of 
variables corresponds to arcs from V + to V~ . 





(i) (ii) 

Figure 5: Cutsets corresponding to binomials (i) and (ii). 



Lemma 5.9 g\l G Q for some p and x T± can be reduced by flvU? *.e. £/ie initial term of g,2 corre- 
sponds to the set of arcs {(k,i r + 1) : k < i r }- 

Proof of Lemma \5. The case of r = 1 is trivial. 

t) 



We suppose r > 1 and x Tl cannot be reduced by any gr\ . Then x Tl can be reduced by some 



which is an element of and x Tl can be also reduced by 9$ t \ (otherwise, g^ t \ is reduced by 

(2) 

and 9(J t ) cannot be an element in Q). 



Suppose that x 1 1 can be reduced by 9^ t \ with t = 1. Let m\ be the monomial obtained by 

reducing x Tl by 9^ t \ , then m\ can be reduced to the monomial 7772 by gjjj (the initial term of g^l is 
a product of variables corresponds to arcs from V ~ to V + by assumption) . 



reduce by 9q\\ 


reduce by g^l 


divided variables 


multiplied variables 


divided variables 


multiplied variables 


{x k ,i r +i : k < i r }, 

{ x k,i qW +l '■ 

k < i q (x), k^i r + l} 


{xir+i,l : 

l>i r + 2, l^i q(1) + l}, 

{x iq(1)+ i,i : I > i q{ i) +2} 


I > i r + 2} 


{xk,i r +i ■ k <i r } 



Table 4: Divided and multiplied variables while reducing by <?n \) and g^ . 

For a binomial fu £ I(i -m t ) which corresponds to the cutset D for (V^,V^) such that V£ = 
{i q (i) + 1} and = V\ , in^(fo) corresponds to arcs from Vp to Vp (otherwise, x T "\M can be 
reduced by fr>, which contradicts the assumption of the induction). Then 777-2 can be reduced by fr>, 
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and the resulting monomial is x Tl (see Table 3), which contradicts to the definition of a term order 
by b. 

Similarly, in the case that x T1 can be reduced by g$ t \ for some t > 1, using fp G 1^ _ m t) 
which corresponds to the cutset D for (Vp, V£) such that V£ = {i q (\) + 1, i g (2) + !)•••) ^q(t) + 1}, and 
= V \ Vp, we can show a contradiction. Thus there exists some p such that gty G 
If a; Tl cannot be reduced by g9^l, i.e. the initial term of gty corresponds to the set of arcs 

{(i r + 1,1) : I > i r + 2}, then can be reduced by gQ, which contradicts that g^ is an element of 
reduced Grobner basis Q. Thus the second statement follows. □ 

We show that if x Tl G in-Al^ -m t ))i then x T2 cannot be reduced by any binomial in Q. If x T2 
can be reduced by some g G G, then such g is one of the following form. 

(i) the binomial g^y and its initial term is Xi r ^ r+ \, 

(ii) any binomial which corresponds to the cutset for (V + , V~) such that i r + 1 G V + and 1,2,..., i r , 

i r + 2 G V~, and its initial term is a product of variables correspond to arcs from V + to V~ , or 

(iii) (The case of r > 1) g^ t y and its initial term is a product of variables correspond to arcs from 
V+ to V~. 




Figure 6: Cutsets corresponding to binomials (i), (ii) and (iii). 



If the case (i) occurs, the initial term of g^K is Xj ri i r +i, which contradicts that i r ^ M(0). On 

the other hand, a binomial of type (ii) can be reduced by g^ by the above lemma, and cannot be 
contained in Q. 

Let us consider that the case (iii) occurs. If x TZ can be reduced by gff ^ with t = 1. Then the 

monomial to which x Tl are reduced by g^\^ can be reduced by a binomial fo G /(/ -m t )i f° r the 
cutset D which corresponds to (V^, Vp) where V£ = {1, 2, . . . , i r — 1, i r + 1} and Vp = V \ V£ , to 
some monomial m (the initial term of fu is a product of variables corresponds to arcs from v£ to Vp 
since i r <£ M(0)). 

For a binomial G _m t ) which corresponds to the cutset D' for (V^,, V^",) such that Vp, = 
{i q (i) + 1} and V£, = V \ V^,, m^/o') corresponds to arcs from V^, to V^, (otherwise, ar 7 ^^} can 
be reduced by fo', which contradicts the assumption of the induction). Then m can be reduced by 
fo' , and the resulting monomial is x T2 (see Table 4) , which contradicts to the definition of a term 
order by b. 
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reduce by g { ^ A) 


reduce by fp 


divided variables 


multiplied variables 


divided variables 


multiplied variable 


X ir,ir+ll X ir,iq(l)+ll 
2 'jr+2,i q ( 1 )+l) 

x i r +3,i q ^+lt 

• • • ' X *g(l) A7(l)+1 


{ x k,i r '■ k < i r — 1}, 
{ x k,l '• k < v + 1, k ^ i r , 
l>i r + 2,l^ + 1}, 

favD+i,! : * > Mi) + 2 } 


{ x k,i r '• k < v — 1}, 
{x M : fc < i r + 1, 
k ^ i r ,l > i r + 2} 


X ir,ir + l 



Table 5: Divided and multiplied variables while reducing by g^ ^ and fo- 

Similarly, in the case that x T2 can be reduced by gfft) f° r some t > 1, using the same fo and 
fo' € /(J _m t ) which corresponds to the cutset D' for (V^,, V^,) such that V^, = {i q (i) + Mg(2) + 
1, . . . , i q ^ + 1}, and V^, = V \ V D ,, we can show the contradiction. 

Therefore, x T2 ^ in^(I^ _m t ))i and T 2 is a desired spanning tree T a . □ 
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